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摘要 : 非 线 性 介质 是 指 其 电导 率 和 (或 ) 介 电 常数 随 电场 强度 变化 而 变化 的 电 介 
质 。 非 均匀 电场 中 非 线性 介质 由 于 介质 的 介 电 性 能 参数 在 空间 呈 梯 度 变化 ， 在 电场 
作用 下 将 发 生 空 间 电荷 极 化 ， 而 极 化 形成 的 空间 电荷 又 影响 电场 的 分 布 ， 探 明 这 种 
空间 电荷 极 化 的 动态 过 程 有 助 于 高 压 直 流 绝缘 结构 的 设计 和 故障 分 析 。 为 此 ， 利 用 
COMSOL Multiphysics 有 限 元 软件 中 的 瞬 态 求解 器 仿真 研究 了 同 轴 电 极 结构 下 具有 场 
致 增强 型 和 温度 增强 型 电导 的 非 线性 电介质 在 阶 跃 电压 下 的 极 化 和 退 极 化 过 程 。 研 究 
结果 发 现 : 极 化 建立 过 程 中 ， 内 电极 附近 空间 电荷 密度 随时 间 变 化 会 出 现 “ 过 冲 ” 现 
象 ， 退 极 化 过 程 中 ， 内 电极 附近 空间 电荷 会 出 现 “ 反 极 性 过 冲 ” 现 象 ;“ 过 冲 ” 现 象 与 
“ 反 极 性 过 冲 ” 现 象 相 互 呼应 ， 均 源 自 动态 松弛 时 间 的 时 空 分 布 ， 温度 梯度 与 场 强 对 空 
间 电荷 “过 冲 ” 行 为 的 影响 规律 不 同 。 

天 键 词 : 非 线 性 介质 ”温度 梯度 ”空间 电 葡 ” 极 化 ” 退 极 化 
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Abstract: The nonlinear dielectric is a sort of dielectrics that their conductivity 
and (or) permittivity will change with the electric field. Under inhomogeneous field. 
The nonlinear dielectric will polarize for its conductivity and (or) permittivity that have 
strong dependence on temperature and electric field. The space charge accumulated in 
polarization process will affect the electric field pattern. To ascertain the characteristics 
of this space charge polarization is of a great significance to help to design the insulation 
structure and to diagnose the insulation breakdown. In this paper, a coaxial model of high 
voltage direct current (HVDC) cable insulation was built and the finite-element analysis 
software (COMSOL Multi-physics) was employed to simulate space charge dynamic 
process caused by the nonlinear conductivity of insulation with different temperature 
gradients and step voltage amplitudes. The results show that a space charge density 
overshoot phenomenon occurs during polarization and an opposite polarity charge density 


overshoot occurs during depolarization near the inner electrode. The polarity charge 
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density overshoot and opposite polarity charge density overshoot are both induced by the distribution of dynamic 


relaxation time in dielectrics. The phenomenon of overshoot during polarization and depolariza-tion echoes 


each other. The temperature gradient and the electric field have different influence on space charge overshoot 


phenomenon. 


Keywords: Non-linear dielectrics, temperature gradient, space charge, polarization, depolarization 


1 引言 


随 着 柔性 输电 技术 的 发 展 ， 聚 合 物 挤 出 绝缘 高 
压 直流 电缆 具有 广泛 的 应 用 前 景 "…"。 众 所 周知 ， 直 
流 稳 态 电场 分 布 取决 于 绝缘 材料 的 电导 率 ， 因 此 在 
存在 内 高 外 低 的 温度 梯度 条 件 下 ， 聚 合 物 绝缘 高 压 
直流 电缆 中 由 于 绝缘 材料 电导 非 线性 因素 会 引起 电 
缆 绝 缘 中 出 现 “电场 分 布 翻转 ”现象 ， 即 电缆 线 世 
表面 电场 强度 低 于 绝缘 外 表面 处 电场 强度 的 现象 ”|。 
由 于 最 大 电场 是 绝缘 结构 设计 的 主要 依据 ， 工 程 界 
往往 只 关注 电场 分 布 问题 。 然 而 “电场 分 布 翻转 ” 
的 本 质 是 绝缘 介质 在 直流 电场 作用 下 慢 极 化 形成 的 
空间 电荷 所 引起 的 。 

长 期 以 来 ， 有 关 空 间 电 荷 的 产生 、 聚 积 和 耗 散 
机 理 以 及 空间 电荷 测试 和 空间 电荷 抑制 方面 的 研究 
一 直 是 电气 绝缘 领域 的 研究 热点 “。 绝 缘 介质 中 
空间 电荷 往往 来 源 于 绝缘 介质 中 杂质 的 热电 离 、 极 
性 分 子 的 极 化 、 电 极 的 注入 、 陷 阱 、 介 电 常 数 或 电 
导 在 绝缘 介质 中 的 不 均匀 分 布 以 及 界面 的 影响 下。 
为 此 可 将 空间 电 共 分 为 两 类 ， 一 类 是 产生 在 强 电场 
下 取决 于 绝缘 材料 微观 特性 (如 陷阱 等 因素 ) 的 空 
间 电 和 荷 ， 另 一 类 是 由 材料 宏观 介 电 参数 在 空间 变化 
而 形成 的 慢 极 化 引起 的 空间 电荷 。 针 对 前 者 ， 相 关 
研究 报道 较 多 "J， 但 引起 高 温度 梯度 下 高 压 直 流 
电缆 电场 分 布 翻转 的 正 是 后 一 类 空间 电荷 ， 它 与 绝 
缘 材 料 电导 率 非 线性 属性 、 结 构 参 数 、 温 度 梯度 和 
外 施 电压 等 多 种 因素 有 关 ， 且 发 生 在 任何 电场 下 。 
由 于 空间 电荷 测试 技术 无 法 区 分 这 两 类 空间 电荷 ， 
通常 认为 后 一 种 空间 电荷 数值 上 也 相对 较 小 ， 
向 常 被 研究 者 忽略 ， 至 今 相 关 研 究 报道 也 较 少 。 

仿照 高 压 直流 电缆 绝缘 结构 ， 本 文 建立 二 维 轴 
对 称 同 轴 电 极 结构 模型 ， 基 于 COMSOL 有 限 元 分 
析 软 件 ， 利 用 数值 法 求解 绝缘 介质 中 电场 时 空 分 布 ， 
进而 得 到 极 化 、 退 极 化 过 程 中 空间 电荷 时 空 分 布 的 
变化 规律 。 研 究 结 果 有 助 于 对 复杂 绝缘 结构 中 非 线 
性 绝缘 介质 空间 电荷 极 化 特性 的 理解 。 


2 同 轴 电极 结构 下 非 线性 介质 空间 电荷 
极 化 仿真 模型 


2.1 仿真 模型 

参考 320kV，500MW 直流 电缆 典型 结构 3， 
对 其 进行 简化 后 建立 如 图 1 所 示 的 二 维 轴 对 称 同 轴 
电极 结构 模型 ， 模 型 中 内 电极 半径 为 26mm， 非 线 
性 绝缘 介质 厚度 为 24mm， 外 电极 半径 为 52mm。 


多 


电极 ( 铀 ) 


非 线性 介质 
(计算 区 域 ) 
内 电极 ( 铜 ) 


图 1 同 轴 电 极 结构 下 非 线性 介质 仿真 模型 
Fig.l Simulation model of nonlinear dielectrics in 


coaxial electrode 


2.2 材料 属性 、 边 界 条 件 与 加 载 方式 
2.2.1 材料 属性 
模型 求解 时 各 部 分 材料 属性 见 表 1。 
表 1 模型 各 部 分 材料 属性 
Tab.1 Material properties in the model 


内 、 外 电极 ( 铜 ) 介质 
5.998 x 107 y(E, T) 


HL SERS y/(S m’) 


相对 介 电 常数 E, 108 2.3 


表 1 PETERE E, T) 由 文献 24] 给 出 ， 即 


sinh (BJE D 


y(E,T) = Aexp 
| E | 


(1) 


_ PA 
k 


b 


式 中 , /为 电导 率 (Sm '); 4 为 与 材料 有 关 的 常 
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数 (V: (Q - m)') ; Pe 为 活化 能 (eV); q AITE E T 
量 ; 有 6 为 玻 尔 兹 曼 常 数 (IK), 7 为 材料 绝对 温度 
(K); B 为 电导 率 对 电场 的 依赖 系数 (m-V'); EA 
电场 强度 (V .m )。 由 实测 电导 率 数据 按照 式 (6) 
进行 数据 拟 合 ， 得 到 4 =3.59 x 10'V (Q: my", g.= 
0.96eV, B= 1.14x10’m-V", 
2.2.2 温度 场 的 确定 
模型 中 温度 边界 条 件 的 设 定 见 表 2。 


#2 温度 边界 条 件 


Tab.2 Boundary conditions of temperature 


边界 URE /K 
AHL FT, Toy t AT 
外 电极 边界 Tou 313 


通过 改变 表 2 中 AT 实现 介质 中 不 同 的 温度 梯 
EDT, P r 处 介质 的 温度 为 


Inr-InR 
TO =T, -—— hn 


六 (2) 
m InR -InR， 


AP, Ra Roa 分 别 为 内 、 外 电极 半径 , 为 介质 中 
任意 一 点 处 半径 。 
2.2.3 电场 边界 条 件 的 确定 

外 电极 接地 ， 内 电极 施加 由 式 (3) 确定 的 激励 
电压 。 


0 t<0s 
U(t)=1 U, Os<t<50000s (3) 
0 50000s <t 


研究 电压 幅 值 对 非 线 性 介质 空间 电荷 极 化 的 
影响 时 ， 分 别 设 置 6 为 正极 性 100kV、200kV.、 
300kV, 400kV 和 500kV。 

2.3 数值 求解 方法 

采用 软件 COMSOL Multiphysics 有 限 元 软件 中 
的 瞬 态 求解 器 对 模型 求解 ， 其 本 质 是 对 泊 松 方程 
进行 数值 求解 过 程 。 数 值 求解 过 程 中 ， 空 间 网 格 
剖 分 单元 大 小 和 时 间 步 长 是 影响 求解 精度 的 关键 
因素 。 对 计算 区 域 进 行 三 角 网 格 剖 分 时 ， 剖 分 单元 
最 大 尺寸 设 定 值 为 0.2mm。 为 了 充分 保证 极 化 和 退 
极 化 过 程 能 达到 稳 态 ， 计 算 时 间 为 150 0008， 其 中 
前 50 000s 为 极 化 过 程 ， 后 100 000s 为 退 极 化 过 程 。 
为 了 保证 计算 精度 和 缩短 计算 时 间 ， 采 取 分 时 段 设 
定 计算 步 长 ， 有 具体 设 置 见 表 3。 
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表 3 仿真 计算 步 长 


Tab.3 Simulation time step 


物理 过 程 时 间 段 4 步 长 /s 
0 ~ 1(50 000 ~ 50 001) 0.01 
1 ~ 10(50 001 ~ 50 010) 0.1 
极 化 过 程 10 ~ 100(50 010 ~ 50 100) 1 
(括号 内 为 退 极 化 过 程 ) ”100 ~ 1 000(50 100 ~ 51 000) 10 


1000 ~ 10 000(51 000 ~ 60 000) 100 


10 000 ~ 50 000(60 000 ~ 150000) 1 000 


2.4 数值 求解 结果 

COMSOL Multiphysics 有 限 元 软件 计算 结果 可 
直接 提供 电势 gq、 电场 强度 EE、 电 位 移 D 和 电 答 密 
Ep 的 时 空 分 布 ， 但 直接 提取 的 空间 电荷 密度 p 的 
时 空 分 布 结果 存在 模型 误差 (由 剖 分 网 格 精细 程度 
决定 )， 电 极 边界 处 空间 电荷 密度 明显 和 物理 事实 不 
符 。 为 此 ， 本 研究 采取 由 电位 移 DD 的 整体 时 空 分 布 
推算 出 电极 附近 受 模 型 误差 影响 处 的 电位 移 人 ， 进 
一 步 计 算出 空间 电荷 密度 p 的 时 空 分 布 ， 具体 关系 
式 为 


(4) 


式 中 ,，D; 和 Di 分 别 为 半径 7 与 ri 处 对 应 的 电 
位 移 。 


3 空间 电荷 极 化 建立 过 程 


3.1 空间 电荷 极 化 建立 过 程 的 温度 梯度 效应 

为 了 研究 温度 梯度 对 极 化 建立 过 程 中 空间 电荷 
的 影响 规律 ， 在 内 电极 施加 电压 幅 值 U = 300K 的 
正极 性 阶 跃 电压 ， 分 别 在 温度 梯度 为 AT= OK、AT7= 
10K, AT=20K, AT=30K, AT=40K fl AT=50K 
条 件 下 进行 了 电场 数值 仿真 ， 仿 真得 到 电场 强度 时 
空 分 布 云图 如 图 2 所 示 。 

从 图 2 可 以 看 出 ， 随 着 温度 梯度 的 增加 ， 稳 态 
电场 分 布 逐 渐 出 现 翻 转 。 由 于 边界 条 件 没 变 ， 介 质 
内 部 电场 发 生 翻 转 的 本 质 必 然 是 空间 电荷 的 积累 。 
通过 计算 得 到 与 图 2 各 个 温度 梯度 下 电场 分 布 相对 
应 的 空间 电荷 密度 分 布 云图 ， 如 图 3 所 示 。 

从 图 3 可 以 看 出 ， 在 正极 性 阶 跃 电 压 作 用 下 ， 
极 化 过 程 介质 中 积累 的 空间 电荷 均 为 正极 性 空间 电 
Pi FAEH (阳极) 来 说 为 同 极 性 空间 电荷 ， 
将 会 削弱 内 电极 附近 电场 ， 对 于 外 电极 〈 阴 极 ) 而 
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2 不 同 温度 梯度 下 电场 强度 时 空 分 布 云图 
Fig.2 Temporal and spatial distribution of electric field with 


different temperature gradients 


言 为 异 极 性 空间 电荷 ， 将 会 增强 外 电极 附近 电场 ， 
使 得 电场 逐渐 发 生 如 图 2 所 示 的 翻转 变化 。 在 所 研 
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(f) AT=50K 
图 3 不 同 温度 梯度 下 介质 中 的 空间 电荷 时 空 分 布 云图 


Fig.3 Temporal and spatial distribution of space charge with 


different temperature gradients 
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究 电压 幅 值 下 ， 稳 态 径 向 空间 电荷 密度 分 布 形 式 随 
着 温度 梯度 的 增加 逐渐 由 内 高 外 低 转变 为 内 低 外 高 。 
在 相近 的 温度 和 场 强 条 件 下 ， 空 间 电 珍 密 度数 值 与 
文献 [13]、 文 献 [22]、 文 献 [25] 和 文献 [26] 实测 数 
值 在 同一 数量 级 。 

介质 内 稳 态 平均 空间 电荷 密度 可 以 反映 其 稳 态 
时 积累 的 空间 电荷 总 量 ， 介 质 的 稳 态 平均 空间 电荷 
密度 与 温度 梯度 关系 如 图 4 所 示 。 从 图 4 可知， 在 
确定 电压 下 钦 度 梯度 越 大 ， 平 均 空 间 电荷 密度 越 大 。 


0.030 | 
0.025 | 
0.020 | 
0.015 + 
0.010 F 
0.005 上 
0.000 


s fl a 7228 BE/(C-m*) 


Be 
S 


0 10 20 30 40 50 
温度 梯度 到 

4 ” 稳 态 平均 空间 电荷 密度 与 温度 梯度 关系 

Fig.4 Steady average space charge density vs. 


temperature gradient 


在 图 3 中 ，A7T=0K 时 径 向 电荷 密度 随时 间 增 加 
基本 呈现 单调 变化 ， 而 存在 温度 梯度 时 内 电极 附近 
电荷 密度 随时 间 变 化 出 现 “ 过 冲 ” 现 象 ， 即 空间 电 
荷 密度 随 极 化 时 间 先 增加 后 减 小 ， 最 后 趋 于 稳 态 值 
的 非 单调 现象 。 

为 了 更 清晰 地 体现 不 同 温度 梯度 条 件 下 非 线 性 
介质 极 化 过 程 中 出 现 的 空间 电荷 密度 “过 冲 ” 行 为 
的 动态 过 程 ， 分 别 在 介质 r=26mm (WER), r= 
38mm 和 r=50mm (外 电极 ) 位 置 处 提取 空间 电荷 
密度 ， 得 到 空间 电荷 密度 随 极 化 时 间 变 化 的 曲线 ， 
如 图 5 所 示 。 

空间 电荷 极 化 “过 冲 ” 现 象 很 不 明显 。 随 着 温 
度 梯 度 增 大 ， 空 间 电荷 密度 “过 冲 ” 现 象 逐 渐 凸 显 
出 来 ， 任 意 位 置 空间 电荷 达到 稳 态 的 时 间 缩 短 ， 内 
电极 暂 态 空间 电荷 峰值 增 大 ， 在 介质 中 部 位 置 暂 态 
空间 电荷 过 冲 现象 逐渐 销 失 ;在 外 电极 附近 极 化 过 
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5 不 同 温度 梯度 下 不 同 部 位 介质 中 空间 电荷 极 化 过 程 
Fig.5 Space charge polarization in dielectrics at different 


positions under different temperature gradients 


程 则 呈现 单调 性 。 同 时 ， 内 电极 附近 稳 态 空间 电荷 
密度 随 着 温度 梯度 的 增加 先 增 大 后 减 小 ， 外 电极 附 
近 空 间 电荷 密度 单调 增加 。 

极 化 空间 电荷 的 产生 是 介质 电导 率 在 空间 呈 梯 
度 分 布 的 结果 ， 不 同位 置 电荷 密度 随时 间 的 变化 过 
程 由 式 (5) 决定 


pentran PED TY) (5) 


AP, er) 为 介质 中 载 流 子 的 动态 松弛 (响应 ) 时 
间 ， 由 下 式 表示 


£ 
HW (6) 
AF, 为 介 电 第 数 ， 与 温度 和 电场 无 关 ; yr) 为 
介质 电导 率 的 时 空 分 布 。 
空间 电荷 的 动态 行为 很 大 程度 取决 于 松弛 时 间 
的 时 空 分 布 情况 。 由 计算 得 到 不 同 温度 梯度 下 松弛 
时 间 tr) 的 时 空 分 布 如 图 6 所 示 。 
从 图 6 可 以 看 出 ， 介 质 的 动态 松弛 时 间 在 任意 
时 刻 沿 半径 方向 均 呈 现 内 低 外 高 的 分 布 形 式 ， 施 加 
电压 之 初 ， 动 态 松弛 时 间 的 空间 分 布 梯度 高 于 稳 态 
时 动态 松弛 时 间 空 间 分 布 梯度 。 在 内 电极 施加 正极 
性 电压 时 ， 带 正 电 的 载 流 子 由 内 电极 向 外 电极 迁移 ， 
由 于 内 导体 附近 7 处 松弛 时 间 小 于 r+Ar 处 松弛 时 
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12000 累 反 过 来 降低 了 > 处 的 电场 ， 并 使 电导 率 下 降 ， 载 


eae ee 
6000 达到 稳 态 。 

4000 由 图 6 同样 可 知 ， 随 着 温度 梯度 的 增 大， 动态 
ie 松弛 时 间 的 时 空 梯度 增 大 。 这 将 导致 施加 电压 初期 


内 电极 附近 电荷 的 积累 速率 过 快 ， 进 而 导致 内 电极 
附近 空间 电荷 极 化 过 程 中 表现 出 “过 冲 ” 现 象 ， 且 
温度 梯度 越 大 ,“ 过 冲 ” 现 象 越 明 显 。 介 质 松 弛 时 间 


由 内 向 外 逐渐 增 大 ， 相 应 电荷 积累 速率 逐渐 减 小 ， 
因此 空间 电荷 表现 出 如 图 5 所 示 规 律 。 

为 了 表述 空间 电荷 极 化 暂 态 过 程 的 “过 冲 ” 行 
为 ， 以 AT=30K、U=300kV 时 结果 为 例 ， 提 取 极 


松弛 时 间 /s 
S 
S 


dp 
35° 7901 0'102 0 
ss 30; 25 {9210 110° O eits 


(b) AT= 10K 化 过 程 中 暂 态 最 BK HEL fay 密度 和 稳 态 态 电 Aap TE 度 沿 半径 
ea 方向 的 分 布 ， 如 图 7 所 示 。 
量 a fg 0.024} KESE 一 暂 态 最 大 值 ， 
= 3000- © 0.022} ; 
S p a x 0.020} | : 
40 35 11021010 = 0.018} ;外电 极 | 
Fim 3 2510" io cate $ 0016| : 
m f i 
(c) AT=20K Za 0.014} : 
0.012 i 
12000 25 30 35 40 45 50 
12000 10000 半径 /mm 
= la 8000 7 ， 径 向 稳 态 与 最 大 暂 态 空间 电荷 密度 分 布 
bz a a on Fig.7 Distribution of maximum transient and steady space 
50 45 40 ie 1010710 31010 5999 charge density along the radius 
Fim 30 25 19°10 fio vest 5 
T 由 于 仿真 结果 存在 模型 误差 ， 图 7 中 两 条 曲线 
d) AT=30K 
00。 RA AE EEN, UE OMS 
a 12000 .4 10000 态 值 的 差 值 小 于 稳 AAW 3 定义 为 临 界 点 wr) 以 此 
= 9000 My 8000 认为 两 点 重合 ， 标 记 为 A (A 点 的 确定 标准 并 不 
里 3000~ if 影响 后 续 对 结果 的 分 析 ) ，A 点 到 内 电极 之 间 为 “过 
10°10" f 2000 冲 ” 区 域 ， 同时 以 内 电极 边缘 处 电荷 密度 暂 态 最 大 
30. 10° 加 
tig, S100 so 0 值 与 对 应 位 置 稳 态 值 的 差 与 稳 态 值 之 比 上 表示 最 大 
(e) AT=40K “过 冲 ” 的 程度 。 
a we a 基于 上 述 定义 ， 得 到 确定 的 电压 幅 值 ， 不 同 温 
‘S000 SS 度 梯度 下 ， 极 化 过 程 中 与 空间 电荷 “过 冲 ” 程 度 和 
ao “过 冲 ” 区 域 ， 如 图 8 所 示 。 由 图 8 可 知 ， 在 确定 的 
| + — ` » N SH: e yh? 
$a 电压 幅 值 下 ， 温 度 梯度 的 增加 使 得 空间 电荷 “过 冲 
llm 30251910 as 程度 增 大 ,“ 过 冲 ” 区 域 变 宽 。 
(i) AT=50K 3.2 空间 电荷 极 化 建立 过 程 的 电压 幅 值 效应 
zz] SHEO HA + +o 4: 、 = 图 > R ye 
ee 为 了 研究 电压 幅 值 对 空间 电荷 极 化 过 程 的 影响 ， 
Fig.6 Temporal and spatial distribution of the relaxation 在 固定 温度 梯度 AT= 50K 的 前 提 下 ， 对 幅 值 分 别 为 
time in polarization 100kV, 200kV, 300kV, 400kV 和 500kV 的 正极 性 


间 , 处 电荷 迁移 速率 大 于 r+Ar 处 电荷 迁移 速率 ， 阶 跃 电压 作用 下 的 极 化 过 程 进行 了 仿真 ， 得 到 极 化 
因而 形成 正极 性 空间 电荷 的 积累 ， 同 时 正 电 荷 的 积 。 ”过程 介质 中 典型 电场 时 空 分 布 去 图， 如 图 9 所 示 。 
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(b) 
图 8 极 化 过 程 温度 中 梯度 对 空间 电荷 过 冲程 度 与 
区 域 的 影响 


Fig.8 The degree and the range of space charge overshoot in 


过 冲 区 域 /mm 
ouav BPa o 


polarization related with the temperature gradient 
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(b) U=500kV 
9 不 同 电压 幅 值 下 电场 强度 时 空 分 布 云图 
Fig.9 Temporal and spatial distribution of electric field with 


different voltages 


结合 图 9 与 前 文 图 2f (AT=50K, U=300kV 时 
的 电场 时 空 分 布 云图 ) 可 以 发 现 ， 电 压 幅 值 越 大 ， 
电场 翻转 过 程 越 快 。 

同样 ， 通 过 计算 可 得 到 与 图 9 中 不 同 电压 幅 值 
下 电场 时 空 分 布 相对 应 的 空间 电荷 密度 时 空 分 布 云 
图 ， 如 图 10 所 示 。 
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(b) U=500kV 
10 不 同 电压 幅 值 下 空间 电荷 密度 时 空 分布 云 图 
Fig.10 Temporal and spatial distribution of space charge 


density with different voltages 


急 态 时 介质 的 平均 空间 电荷 密度 与 极 化 电压 关 
系 如 图 11 所 示 。 由 图 11 可 知 ， 介 质 中 空间 电 衔 总 
量 随 极 化 电压 增 大 而 增 大 ， 但 增加 速率 随 极 化 电压 
增 大 而 减 小 。 


0.040 | 
0.035 | 
0.030 | 
0.025 | 
0.020 | 
0.015 + 


25 (i FLL Af 3 RE/(C-m”) 


100 200 300 400 500 
极 化 电压 /kV 


ll 稳 态 平均 空间 电荷 密度 与 极 化 电压 关系 
Fig.ll Steady average space charge density vs. 


polarization voltage 


为 了 更 清晰 地 对 比 不 同 电压 幅 值 下 介质 中 的 空 
间 电 和 荷 极 化 动态 过 程 ， 同 样 做 出 半径 r+=26mm (内 
电极 )、r=38mm 和 x=50mm (外 电极 ) 处 的 空间 
电荷 密度 随 极 化 时 间 变 化 的 曲线 ， 如 图 12 所 示 。 

由 图 12 可 知 ， 随 电压 幅 值 增 大 ， 任 何 位 置 空间 


g =y: 
5 0.04 =(= 
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图 12 不 同 电压 幅 值 下 不 同 部 位 介质 空间 电荷 极 化 过 程 
Fig.12 Space charge polarization in dielectrics at different 


positions under different voltage amplitudes 


电荷 达到 稳 态 的 时 间 缩 短 ， 内 电极 暂 态 空间 电荷 峰 
值 增 大 ;在 介质 中 部 暂 态 空间 电荷 过 冲 现 象 逐 渐 消 
Fes 外 电极 附近 空间 电荷 建立 过 程 无 过 冲 现象 ， 且 
稳 态 电荷 密度 随 极 化 电压 增加 而 增加 。 

在 确定 的 温度 梯度 条 件 下 ， 极 化 过 程 中 电压 幅 
值 对 空间 电荷 “过 冲 ” 程 度 和 “过 冲 ” 区 域 的 影响 
如 图 13 所 示 。 

由 图 13b 可 知 ， 当 A7T=0K 时 ,在 电压 大 于 
等 于 300kV 时 才 出 现 明显 “过 冲 ” 现 象 ， 而 后 随 
着 电压 幅 值 的 增加 “过 冲 ” 程 度 与 区 域 增 大 。 当 
AT > 10OK 时 ， 电 压 幅 值 的 增 大 使 得 “过 冲 ” 程 度 
与 区 域 减 小 。 由 此 可 见 ， 电 压 幅 值 对 极 化 “过 冲 ” 
现象 影响 规律 相反 。 由 式 (1) 电导 率 公式 可 知 ， 温 
度 梯度 越 高 ， 电 导 率 空间 分 布 梯度 越 大 ， 而 电压 幅 
值 的 增加 使 同 轴 结 构 中 场 致 增强 型 非 线 性 绝缘 介质 
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13” 极 化 过 程 中 电压 幅 值 对 空间 电荷 过 冲程 度 与 
区 域 的 影响 


Fig.13 The influence of voltage on the degree and the range 


of space charge overshoot in polarization 


rH EL Sy > is AP), Emh LS 8 8 Td Ah 
梯度 。 


4 空间 电荷 退 极 化 耗 散 过 程 


4.1 空间 电荷 退 极 化 耗 散 过 程 的 温度 梯度 效应 

研究 温度 梯度 对 退 极 化 过 程 空 间 电荷 的 影响 
时 ， 极 化 电压 UV = 300kV， 温 度 梯度 分 别 设置 为 A7 
=0K, AT=10K, AT=20K, AT=30K, AT=40K 和 
A7T=50K， 在 上 述 条 件 下 对 空间 电荷 退 极 化 过 程 进 
行 仿真 ， 得 到 空间 电荷 退 极 化 过 程 电 场 时 空 分 布 如 
图 14 所 示 。 

由 图 14 可 知 ， 退 极 化 过 程 中 ， 内 电极 附近 电场 
与 极 化 时 方向 相反 ， 而 外 电极 附近 电场 与 极 化 时 同 
向 ， 中 间 区 域 存 在 一 个 电场 为 零 的 点 ， 而 零 电 场 出 
现 的 位 置 在 退 极 化 过 程 中 是 变化 的 ， 最 终 稳 态 时 电 
场 处 处 为 零 。 

根据 电场 时 空 分 布 计算 得 到 不 同 温度 梯度 下 退 
极 化 过 程 空间 电荷 密度 时 空 分 布 ， 如 图 15 所 示 。 

由 图 15 可 知 ， 在 温度 梯度 A7T=0K 时 ,介质 
中 空间 电荷 密度 随时 间 单 调 减 小 ， 这 与 该 温度 梯 
度 下 极 化 过 程 没 有 出 现 “ 过 冲 ” 现 象 相 吻 合 ， 当 
AT > 10K 时 ， 随 着 温度 梯度 的 增加 ， 内 电极 附近 
空间 电荷 出 现 “ 反 极 性 过 冲 ” 现 象 ， 即 内 电极 附近 
空间 电 符 不 是 随时 间 单 调 减 小 到 零 ， 而 是 出 现 反 极 
性 的 峰值 后 再 衰减 到 零 。 退 极 化 过 程 的 “ 反 极 性 过 
冲 ” 与 极 化 过 程 的 “过 冲 ” 相 互 呼应 ,“ 反 极 性 过 
冲 ” 仍 是 由 动态 松弛 时 间 的 时 空 分 布 所 决定 的 。 由 
于 退 极 化 过 程 内 电极 附近 出 现 反 向 电场 ， 故 出 现 
“ 反 极 性 过 冲 ”。 
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图 14 不 同 温度 梯度 退 极 化 过 程 介质 中 电场 强度 时 空 分 布 Fig.15 Temporal and spatial distribution of space charge in 


Fig.14 Temporal and spatial distribution of electric field in the non-linear dielectric during depolarization with different 


the non-linear dielectric during depolarization with different temperature gradients 


Paranal 为 了 表征 退 极 化 “ 反 极 性 过 冲 ” 程 度 ， 提 取 退 
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极 化 起 始 时 电荷 密度 和 和 暂 态 最 大 电 奏 密度 分 布 ， 结 
采 如 图 16 所 示 。 
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zo OOF 内 电极 外 电极 
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-0.01 ¢ 


z% 


半径 /mm 
16 ， 径 向 初始 与 最 大 暂 态 空间 电荷 密度 分 布 
Fig.16 Distribution of initial and maximum transient space 


charge density along the radius 


图 16 中 ,A 点 为 出 现 “ 反 极 性 过 冲 ” 的 临界 
点 ，A 点 与 内 电极 间 定 义 为 “ 反 极 性 过 冲 ” 区 域 ， 
以 内 电极 边缘 处 负电 荷 密度 峰值 绝对 值 与 对 应 位 置 
初始 值 之 比 左 定义 最 大 “ 反 极 性 过 冲 ” 程 度 。 由 此 
得 到 确定 的 电压 幅 值 下 温度 梯度 对 退 极 化 过 程 空间 
电荷 “ 反 极 性 过 冲 ” 程 度 和 区 域 的 影响 规律 ， 如 图 
17 所 示 。 
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(b) 
图 17 退 极 化 过 程 中 温度 梯度 对 空间 电荷 过 冲程 度 与 
区 域 的 影响 
Fig.17 The influence of temperature gradient on the degree 


and the range of space charge overshoot in depolarization 
由 图 17 可 知 ， 退 极 化 过 程 中 温度 梯度 对 空间 电 


和 荷 “ 反 极 性 过 冲 ” 程 度 和 区 域 的 影响 规律 与 极 化 过 
程 相同 ， 即 “ 反 极 性 过 冲 ” 程 度 与 区 域 均 随 着 温度 
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梯度 的 增加 而 增 大 。 

4.2 空间 电荷 退 极 化 耗 散 过 程 的 极 化 电压 幅 值 效应 
固定 温度 梯度 A7= 50K， 在 极 化 电压 幅 值 分 别 

为 100kV 和 SOOKY 条 件 下 ， 仿 真得 到 此 时 典型 电场 

分 布 云图 与 对 应 的 空间 电荷 密度 分 布 云图 ， 分 别 如 

图 18 和 图 19 所 示 。 
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场 强 /kV.mmD) 


pe T Ak 
(b) U=500kV 

18 不 同 极 化 电压 退 极 化 过 程 电场 强度 时 空 分布 云 图 

Fig.18 Temporal and spatial distribution of electric field 


with different voltages during depolarization 
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(b) U=500kV 
19 不 同 极 化 电压 退 极 化 过 程 空间 电荷 密度 时 空 分 布 


Fig.19 Temporal and spatial distribution of space charge 


density with different voltages during depolarization 


由 图 18 与 图 19 可 知 ， 极 化 电压 增 大 ， 初 始 退 
极 化 电荷 密度 和 退 极 化 过 程 的 电场 随 之 增 大 。 

退 极 化 过 程 中 ， 极 化 电压 幅 值 对 空间 电荷 “ 反 
极 性 过 冲 ” 的 影响 如 图 20 所 示 。 

由 图 20a 可 知 ， 在 空间 电荷 退 极 化 耗 散 过 程 中 ， 
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(b) 
20” 退 极 化 过 程 中 极 化 电压 幅 值 对 空间 电荷 过 冲程 度 
与 区 域 的 影响 


Fig.20 The influence of polarization voltage on the degree 


and the range of space charge overshoot in 


depolarization process 


“ 反 极 性 过 冲 ” 程 度 随 着 极 化 电压 幅 值 的 增 大 而 减 
小 ， 这 与 极 化 过 程 规律 一 致 。 We i 
变 对 空间 电荷 “过 冲 ” 区 域 基本 没有 影响 ， 这 
与 极 化 过 程 不 同 。 


5 结论 


同 轴 结 构 非 线性 介质 中 存在 温度 梯度 时 ， 在 外 
加 电场 的 作用 下 会 发 生 由 电导 率 梯度 引起 的 空间 电 
荷 极 化 。 针 对 这 种 空间 电荷 极 化 与 退 极 化 行为 ， 本 
文 研究 了 温度 梯度 与 极 化 电压 幅 值 对 其 的 影响 ， 得 
出 如 下 结论 : 

(1) 当 极 化 电压 一 定 ， 随 温度 梯度 的 增 大 ， 极 
化 稳 态 空间 电荷 密度 由 内 高 外 低 分 布 逐渐 过 滤 到 内 
低 外 高 分 布 , 温度 梯度 越 大 ， 介 质 中 平均 空间 电荷 
密度 越 大 ， 极 化 到 稳 态 所 需 时 间 越 短 。 

(2) 空间 电荷 极 化 过 程 中 ， 内 电极 附近 空间 电 
茶会 出 现 “ 过 冲 ” 现 象 ， 退 极 化 过 程 中 出 现 “ 反 极 
过 冲 ” 现 象 。“ 过 冲 ” 与 “ 反 极 性 过 冲 ” 相 互 呼 
应 ， 均 源 自 动态 松弛 时 间 的 时 空 梯度 分 布 。 

(3) 温度 梯度 与 场 强 对 空间 电荷 “过 
“ 反 极 性 过 冲 ” 行 为 的 影响 效果 相反 。 


E 2 


ap” Fil 


[2] 


[3] 


[6] 


[7] 


[8] 


[9] 


ChinaXiv 合 作 期 


参考 文献 


温 加 良 ， 吴 锐 ， 彭 畅 ， 等 .直流 电网 在 中 国 的 应 用 
前 景 分 析 四 . 中 国电 机 工程 学 报 ，2012，32(13): 
7-12. 


Wen Jialiang, Wu Rui, Peng Chang, et al. Analysis 
of DC grid prospects in China[J]. Proceedings of the 
CSEE, 2012, 32(13): 7-12. 

fal 4 Pk, WR, AE, S. 挤 压 型 高 压 直 流 电缆 
研究 进展 及 关键 技术 述评 四 .高 电压 技术 ，2015， 
41(5): 1417-1429. 


He Jinliang, Dang Bin, Zhou Yao, et al. Reviews on 
research progress and key technology in extruded 
cables for HVDC transmission[J]. High Voltage 
Engeering, 2015, 41(5): 1417-1429. 
John C Fothergill. The coming of age of HVDC 
extruded power cables[C]. Electrical Insulation 
Conference. Philadelphia, USA, 2014: 124-137. 
Mazzanti G, Marzinotto M. Extruded cables for high 
voltage current transmission: advances in research 
and development[M]. Manhattan, USA: Wiley-IEEE 
Press, 2013: 1-18. 
Rongsheng Liu. Long-distance DC electrical power 
transmission[J]. IEEE Transactions on Dielectrics 
and Electrical Insulation, 2013, 20(5): 37-46. 
Ohki Y. Development of XLPE-insulated cable for 
high voltage DC submarine transmission line[J]. 
IEEE Electrical Insulation Magazine, 2013, 29(4): 
65-67. 
Reddy Ch Chakradhar, Ramu T S. On the 
computation of electric field and temperature 
distribution in HVDC cable insulation[J]. IEEE 
Transactions on Dielectrics and Electrical Insulation, 
2006, 13(6): 1236-1244. 
李 忠 华 ， 刘 乐 乐 ， 郑 欢 ， 等 . HVDC 电缆 电场 分 
影响 因素 的 仿真 研究 四. 中 国电 机 工程 学 报 ， 
2016，36(9): 2563-2571. 
Li Zhonghua, Liu Lele, Zheng Huan, et al. 


Simulation on the influence factors of electric field 
distribution in HVDC cable[J]. Proceedings of the 
CSEE, 2016, 36(9): 2563-2571. 

Mahdavi S, Zhang Y, Alquie C, et al. Determination 
of space charge distribution in polyethylene 


samples submitted to 120kV DC voltage[J]. IEEE 


H 


201903.00132v1 


chinaXiv 


[10] 


[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


[17] 


ChinaXiv 合 作 期 刊 


孙 云 龙 等 ” 同 轴 电 极 结构 中 非 线性 介质 空间 电荷 极 化 与 退 极 化 暂 态 过 程 仿真 研究 2018 年 11 月 


Transactions on Electrical Insulation, 1991, 26(1): 
57-62. 

Fabiani D, Montanari G C, Laurent C, et al. 
Polymeric HVDC cable design and space charge 
accumulation. part 1: insulation/semicon interface[J]. 
IEEE Electrical Insulation Magazine, 2007, 23(6): 
11-19. 

Castellon J, Notingher P, Agnel S, et al. Electric field 
and space charge measurements in thick power cable 
insulation[J]. IEEE Electrical Insulation Magazine, 
2009, 25(3): 30-42. 

Zhang Yibo, Christen Thomas, Meng Xing, et al. 
Research progress on space charge characteristics 
in polymeric insulation[J]. Journal of Applied 
Dielectrics, 2016, 6(1): 1-13. 

UOR, ME, EER. 直流 交 联 聚 乙 烯 绝 缘 中 空 
间 电 荷 的 形成 机 理 [J]. 绝缘 材料 ，2010，43(6) 
39-43. 


Qiu Bin, He Jun, Tu Demin. Formation mechanism 
of space charge in DC XLPE insulation[J]. Insulation 
Materials, 2010, 43(6): 39-43. 
GH, RT, Be, $. 
质 的 界面 性 能 研究 进展 [J]. 
38(9): 2455-2464. 


聚合 物 纳米 复合 电 介 
高 电压 技术 ，2012， 


Luo Yang, Wu Guangning, Peng Jia, et al. Research 
progress on interface properties of polymer 
nanodielectrics[J]. High Voltage Engineering, 2012, 
38(9): 2455-2464. 

KEX, BE, ERE, T. ROMS RAM 
空间 电荷 注入 特性 的 比较 [J]. 高 电压 技术 ，2014， 
40(9): 2613-2618 . 

Zhang Yewen, Zhao Hui, Li Zongze, et al. 
Comparison of space charge injection behaviors 
between polythene and polypropylene[J]. High 
Voltage Engineering, 2014, 40(9): 2613-2618. 

周 党 ， 伍 能 成 ， 廖 瑞金， 等. 不 同 交 联 度 交 联 
聚 乙 烯 的 空间 电荷 特征 [如 . 高 电压 技术 ，2013， 
39(2): 294-301. 


Zhou Quan, Wu Nengcheng, Liao Ruijin, et al. Space 
charge characteristics of cross-linked polyethylene 
with different cross-linking degrees[J]. High Voltage 
Engineering, 2013, 39(2): 294-301. 

RRI HF, JEA, S. 表面 修饰 纳米 SiOy 
XLPE 的 电导 电流 和 空间 电荷 特性 四 .高 电压 技 


[18] 


[19] 


[20] 


[21] 


[22] 


[23] 


术 ，2014，40(10): 3268-3275. 

Wu Zhensheng, Ye Qing, Zhou Yuanxiang, et al. 
Conduction current and space charge characteristics 
of Si02/XLPE nanocomposites with nanoparticles 
surface modification[J]. High Voltage Engineering, 
2014, 40(10): 3268-3275. 

PRR, ZA, RER, T. 交 联 副 产 物 对 交 联 
聚 乙 烯 中 空间 电荷 行为 的 影响 四. 高 电压 技术 ， 
2015, 41(11): 2903-2910. 


Zhong Qiongxia, Lan Li, Wu Jiandong, et al. 
The influence of cross-linked by-products on 
space charge behavior in XLPE[J]. High Voltage 
Engineering, 2015, 41(11): 2903-2910. 

BHG, HEC, Bik, S. teh / ERR Ci 
米 复 合 材料 的 空间 电荷 电导 特性 [J]. 高 电压 技术 ， 
2014，40(9): 2661-2667. 


Yan Zhiyu, Han Baozhong, Zhao Hong, et al. Space 
charge and electrically conductive characteristics 
of CB/XLPE nanocomposites[J]. High Voltage 
Engineering, 2014, 40(9): 2661-2667. 

XNE, HA, Ih, SE. 温度 梯度 影响 高 压 直 流 
电 绕 用 交 联 聚 乙 烯 中 空间 电荷 分 布 的 作用 机 理 [J]. 
高 电压 技术 ，2015，41(8): 2665-2673. 

Liu Tong, Fu Mingli, Hou Shuai, et al. Mechanism 
of space charge distribution in XLPE used in HVDC 


cable under temperature gradient[J]. High Voltage 
Engineering, 2015, 41(8): 2665-2673. 

Fabiani D, Montanari G C, Bodega R, et al. The 
effect of temperature gradient on space charge and 
electric field distribution of HVDC cable models[C]. 
IEEE 8th International Conference on Properties and 
Application of Dielectric Materials, Bali, Indonesia, 
2006: 65-68. 

McAllister I W, Crichton G C, Pedersen A. Space 
charge in DC cables[C]. Conference Record of 
the IEEE International Symposium on Electrical 
Insulation, Montreal, Canada, 1996: 661-665. 
MEH, Ete, HSM, SF. LDPE 纳米 复合 介 
质 的 直流 电导 特性 及 其 对 高 压 直 流 电缆 中 电场 分 
布 的 影响 加. 中 国电 机 工程 学 报 ，2014，34(9): 
1454-1461. 


Yang Jiaming, Wang Xuan, Han Baozhong, et al. DC 
conductivity characteristics of LDPE nanocomposites 


and its effect on electric field distribution in HVDC 


201903.00132v1 


chinaXiv 


[24] 


[25] 


cables[J]. Proceedings of the CSEE, 2014, 34(9): 
1454-1461. 

Steve Boggs, Dwight H Damon, Jesper Hjerrild, et 
al. Effect of insulation properties on the field grading 
of solid dielectric DC cable[J]. IEEE Transaction on 
Power Delivery, 2001, 16(4): 456-461. 

Holboll J T, Henriksen M, Hjerrild J. Space 
charge build-up in XLPE cable with temperature 
gradient[C]. Annual Report Conference on Electrical 
Insulation and Dielectric Phenomena. Victoria, 


Canada, 2000: 157-160. 


[26] 


[27] 


ChinaXiv 合 作 期 刊 


Bambery K R, Fleming R J, Hoboll J T. Space 
charge profiles in low density polyethylene samples 
containing a permittivity/conductivity gradient[J]. 
Journal of Physics D: Applied Physics, 2001, 34(20): 
3071-3077. 

McAllister I W, Crichton G C, Pedersen A. Charge 
accumulation in DC cables: a macroscopic approach 
[C]. Conference Record of the IEEE International 
Symposium on Electrical Insulation, Pittsburgh, 


USA, 1994: 212-216. 


